New Experimental Constraints for the Standard Model from Muon Decay 
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The TWIST Collaboration has completed a new measurement of the energy-angle spectrum of 
positrons from the decay of highly polarized muons. A simultaneous measurement of the muon 
decay parameters p, S, and V^S, tests the Standard Model (SM) in a purely leptonic process and 
provides improved limits for relevant extensions to the SM. Specifically, for the generalized left-right 
symmetric model Ksk/P-OCI < 0.020 and {gL/ga)m2 > 578 GeV/c 2 , both 90% C.L. 

PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 



The maximal parity violation of the Standard Model 
(SM) charged weak interaction is empirically based. 
Many natural SM extensions restore parity conserva- 
tion at a higher mass scale with additional weak cou- 
plings. Muon decay is an excellent laboratory to search 
for these couplings because the purely leptonic process 
can be calculated very precisely within the SM. This Let- 
ter presents a high precision measurement of the energy- 
angle spectrum of the positrons emitted in polarized 
muon decay, which provides new limits for the mass and 
mixing angle of the heavy IT in a class of left-right sym- 
metric (LRS) models [H. 

The most general Lorentz invariant, derivative-free 
muon decay matrix element 0] is described by 10 com- 
plex, model-independent couplings (flftj involving left- 
and right-handed leptons (e, p) in scalar, vector, and ten- 
sor interactions (7). In the SM, g\ L = 1, and the other 
nine constants are zero. When only the positron energy 
and direction are measured, the muon decay spectrum 
can be described by four parameters Q , which are bilin- 
ear combinations of the g]^: the Michel parameter p, as 
well as 5, V^, and 77. The differential decay rate is then 
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The neutrino masses are neglected; E max — 52.828 MeV. 
Radiative corrections Q are not explicitly shown, but 
are significant and must be evaluated within the SM to 
a precision comparable to the experiment. The polar- 
ization of the muon from pion decay begins as and 
may evolve over the 2.2 /zs mean lifetime of the muon to 
become at the time of decay. The SM predictions are 
p = 5 = 3/4, = £ = 1, and 77 = 0. Precision measure- 
ments of these parameters test the SM predictions and 
are sensitive to extensions to the SM. 

Prior to the TRIUMF Weak Interaction Symmetry 
Test (TWIST) experiment, p, S, and VZt, were known 
with uncertainties in the range of (3.5-8.5) x 10 -3 In- 
termediate TWIST results have already reduced those 
uncertainties to (0.7-3.8) x 10" 3 0,0- TWIST has now 
realized its goal of making about an order of magnitude 
improvement in each of the parameters. These final re- 
sults of the experiment supersede those in our previous 
publications. 

TWIST used highly polarized positive muons from pi- 
ons decaying on the surface of a graphite production tar- 
get irradiated with 500 MeV protons. The muons were 
transported by the TRIUMF M13 surface muon chan- 
nel to the entrance of a 2 T superconducting solcnoidal 
magnet and were guided by the field along its symmetry 
axis into the detector array Q. A thin (239 pcm) trigger 
scintillator identified muons entering the detector. The 
muons were ranged to stop predominately in a thin metal 
foil at the center of a symmetric stack of high-precision, 
low-mass planar multiwire chambers. There were 6 pro- 
portional and 22 drift chambers, surrounded by helium 
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gas, on each side of the stopping target. Ionization of 
tracks from decay positrons was sampled by the cham- 
bers and drift times were recorded by TDCs. 

With a muon rate of (2-5) xlO 3 s^ 1 , data sets of 10 9 
events could be obtained within a few days. Data sets 
were taken with two foil stopping targets, Ag (30.9 ^m 
thick) and Al (71.6 pm thick), each >99.999% pure. Sets 
were also taken with deliberately altered conditions to as- 
sist in studies of possible systematic errors. A pair of time 
expansion chambers (TECs) Q was inserted upstream of 
the solenoid to determine the incoming muon beam char- 
acteristics. Because they caused multiple scattering and 
hence muon depolarization, the TECs were removed dur- 
ing most data taking. This phase of the experiment was 
completed in 2007. 

Analyzed events were collected into two-dimensional 
(2D) distributions of positron angle and momentum (or 
energy) whose shape depends on the decay parameters. 
These distributions from data were compared to simi- 
lar ones derived from a GEANT3 simulation [loj . Both 
were subjected to essentially the same analysis, allowing 
biases and inefficiencies to be included in an equivalent 
way to reduce the dependence of the result on the spe- 
cific analysis procedure. This places great importance on 
the accuracy and detail of the simulation, which includes 
not only standard physics processes but also a detailed 
description of the beam, magnetic field, geometry, and 
detector response. Decay parameters were extracted by 
fitting the 2D data distribution to that of a base distribu- 
tion of simulated events, plus simulated distributions cor- 
responding to the first derivatives of the spectrum with 
respect to decay parameters (or combinations thereof), 
yielding fit coefficients Ap, A£, and A£6. The decay pa- 
rameters used in the generation of the simulation were 
hidden, so the analysis was "blind" [f|. 

Fourteen data sets were used to extract p and 5, seven 
with each of the Ag and Al targets. Only nine sets were 
used for "P^£; the other five were acquired to test consis- 
tency and systematic effects with altered beam, magnetic 
field, or muon multiple scattering, where the depolariza- 
tion was not optimally controlled. The residuals for the 
fit of one nominal data set in units of standard deviations 
(a) are shown in Fig. Q] A histogram of the residuals in 
the fiducial region, summed over all sets, has a mean of 
-0.003 ±0.005 and a = 1.002 ±0.004. Also shown is the 
range of (p, cos 9) used to determine the decay parame- 
ters. The fiducial cuts are symmetric for upstream and 
downstream decays and were selected to maximize sensi- 
tivity to the decay parameters while reducing systematic 
uncertainties. For all 14 data sets, there were 11 x 10 9 
events, of which 0.55 x 10 9 passed all event selection cri- 
teria and were within the fiducial region. Simulation sets 
were typically 2.7 times larger than the corresponding 
data set. The consistency of the data sets (statistical un- 
certainties only) was assessed by fits to constant means 
for the values of Ap, A<5, and AV^, which gave reduced 
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FIG. 1. Residuals for the fit of one nominal data set to sim- 
ulation in units of standard deviations. The fiducial region is 
outlined. 



X 2 values of 14.0/13, 17.7/13, and 9.7/8, respectively. 

The procedure of fitting the difference of two spectra 
in terms of derivatives also provides a natural tool for 
the evaluation of systematic uncertainties. The simula- 
tion was validated through comparison with observables 
in the data that do not depend on muon decay parameter 
values, and the resulting uncertainties were factored into 
estimates of the systematic uncertainties. The sensitivi- 
ties were obtained from the effect on the decay parame- 
ters when an identified source of systematic uncertainty 
was changed (often by an exaggerated amount) in one of 
the spectra. This was typically achieved with two sim- 
ulated spectra. The systematic uncertainties are listed, 
along with the statistical errors, in Table fl] 

Notable improvements in the systematic uncertain- 
ties for p and S compared to our intermediate results 
[6j were achieved for positron interactions, chamber re- 
sponse, and momentum calibration. The positron inter- 
actions systematic addresses the possible inaccuracy in 
the simulation of reproducing positron energy loss in the 
stopping target and detector elements, primarily due to 
bremsstrahlung, delta-ray production, and ionization. It 
was better constrained by comparisons of identified in- 
teractions observed in the data and in the simulation. 
Chamber response refers to the conversion of drift times 
to spatial information used in track fitting to evaluate 
the momentum and angle of each decay positron. This 
was improved by more precise monitoring and control 
of atmospheric influences that could change the cham- 
ber cell geometry. In addition, a method was devised 
111 ] to calibrate the chambers' space-time relations, for 
each plane, in both data and simulation, thereby reduc- 
ing reconstruction biases. The maximum positron energy 
provides a calibration feature that was used to reduce the 
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TABLE I. Systematic uncertainties and statistical errors for 
p, S, and 
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Target Independent Systematics 


Momentum calibration 


1.2 


1.2 


1.5 




Chamber response 


1.0 


1.8 


2.3 




Radiative corrections, r\ 


1.3 


0.6 


1.2 




Resolution 


0.6 


0.7 


1.5 




Positron interactions^ 


0.5 


0.2 


0.4 




Others 


0.3 


0.3 


0.4 




Depolarization in fringe field 




+15.8, - 


4.0 


Depolarization in stopping 


target 




3.2 




7r decays in beamline 






1.0 




Uncertainties for Ag 


target 






Bremsstrahlung rate 


1.8 


1.6 


0.5 




Ag thickness/stop position 


3.8 


6.4 


0.6 




Statistical 


1.2 


2.1 


4.2 




Uncertainties for Al target 


Bremsstrahlung rate 


0.7 


0.7 


0.3 




Al thickness/stop position 


0.2 


0.8 


0.8 




Statistical 


1.4 


2.4 


3.9 




Weighted Systematic 


2.3 


2.7 


+16.5 - 


6.3 


Weighted Statistical 


1.2 


2.1 


2.9 




Total Error 


2.6 


3.4 


+16.8 - 


6.9 



a excluding bremsstrahlung 



energy scale uncertainty. Since energy loss varies with the 
track angle linearly in l/(cos 9) due to the planar geom- 
etry of the detector, the region near the kinematic end- 
point of 52.8 MeV/c was matched for data and simulation 
for small bins of cos 9. The data-simulation relative en- 
ergy calibration procedure has undergone improvements 
to become more robust to fitting conditions. 

The asymmetry parameter £ is also subject to uncer- 
tainties from these sources, but they are overshadowed 
by uncertainties unique to depolarization, as shown in 
Table Q] Depolarization in the fringe field and in the 
muon stopping target result in < V* and constitute 
the largest contributions to systematic uncertainties for 
■p^£. These uncertainties were improved considerably for 
this analysis compared to the intermediate result [7j . Im- 
provements in the beam steering reduced the fringe field 
depolarization for a nominal data set to only 2.5 x 1CP 3 . 
The uncertainty depends on the accuracy of simulating 
the muon spin evolution as the beam passes through sig- 
nificant radial field components at the solenoid entrance. 
The essential ingredients are an accurate field map and 
precise knowledge of the position and direction of the 
muons, as provided by the TECs. Depolarization in the 
stopping target from muon spin relaxation is assessed 



from the measured time dependence of the asymmetry. 

After revealing the hidden parameters, the results for 
the three decay parameters are consistent with the SM 
predictions. While the generalized matrix element treat- 
ment of Ref. does not constrain the sign of deviations 
from the SM values for p, 5, and £, the product V^d/p is 
constrained to be < 1.0 and is 1.0 in the SM. This quan- 
tity defines the asymmetry between cos 6* = ±1 at the 
maximum decay positron energy. Our decay parameter 
values combined to give V^S/p = 1.00192 ±^00066 ( the 
errors account for significant correlations) , and the initial 
evaluation of V^S/p showed that the value for the Ag 
data was higher than that for Al by 3.8 a. This apparent 
contradiction initiated an exhaustive reconsideration of 
effects that might have been overlooked in the blind anal- 
ysis. The review showed that effects such as p + -+ e + X° 
decays (where X° is a long-lived unobserved particle), 
an incorrect value of the rj parameter, or plausible er- 
rors in the radiative correction implementation were not 
responsible for the unexpected V^S/p value. 

While no obvious mistake was uncovered in the esti- 
mates of the systematic uncertainties previously consid- 
ered, we found that two corrections had been missed. 
A small correction was added for muon radiative decay 
(<1 x 10~ 4 for the Ag data and negligible for Al). An- 
other correction was made for each set to account for a 
difference between the mean muon stopping position for 
data and simulation. We also concluded that the uncer- 
tainties for the two targets were sufficiently different to 
merit dividing the systematic uncertainties into common 
and target-dependent categories. The target independent 
systematics are unchanged from the blind analysis. Sep- 
arate uncertainties for bremsstrahlung were computed, 
and an additional sensitivity to the muon stopping posi- 
tion in the target was added. 

With these changes the central values of p and S de- 
creased from the blind results by 0.00014 and 0.00023, 
respectively. V^t; is unchanged and its error reduced af- 
ter including information from the measurement of 5 in 
the five sets not used for VZt,- All uncertainties changed 
by < 0.00006. The difference between targets for V^S/ p 
is reduced to -1 a, and V^S/p = 1.00179 to'.ooon- The 
revised results are compared to prior results in Fig. [2] 
The values, including the uncertainties from Table 1, are: 

p = 0.74977 ± 0.00012(stat) ± 0.00023(syst); 
S = 0.75049 ± 0.00021 (stat) ± 0.00027(syst); 
= 1.00084 ± 0.00029(stat) tSlCsyst). 

The decay parameters measured by TWIST contribute 
to a larger set derived from other muon decay observ- 
ables that can be analyzed in terms of the weak couplings 
g]^. A global analysis @, [H| imposes V^S/p < 1.0 and 
yields V^S/p > 0.99909 (90% C.L.), compared to the 
pre-TWIST lower limit P^S/p > 0.99682 (90% C.L.) 
[l4j. The global analysis confirms consistency with the 
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FIG. 2. Summary of published central values and total uncer- 
tainties for p, 8, and VZ£, Hrfi O i along with the results 
of this analysis. Vertical lines represent the SM values. 



model no assumptions are made about the ratio of right- 
to left-handed couplings (gn/gL) or the form of the right- 
handed CKM matrix. In this case, the TWIST result for 
p provides the best limit on the mixing angle between 
the light and heavy mass eigenstates, W\ and W%. Our 
limit is \(g R /g L )C\ < 0.020 (90% C.L.), compared to the 
pre-TWIST limit of \(g R /g L )(\ < 0.066. The lower limit 
on the mass of W2 ((gL/ 'gn)^) has been increased from 
400 GcV/c 2 to 578 GeV/c 2 . Coupled constraints on the 
mass for (<7l/3_r) to 2 an d the mixing angle are shown in 
Fig. [3] where our limits are derived from a correlated 2D 
probability distribution from our measurements. These 
improved constraints will significantly impact predictions 
from the class of LRS models where the neutrinos are 
light compared to the muon mass. 
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of Science, and the U.S. Department of Energy. 
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FIG. 3. Allowed region (90% C.L.) of mixing angle (£) and 
heavy W mass (ni2) for the general LRS model. 



SM, where g\ L is the only non-zero term. The TWIST 
results restrict the upper limits of other terms. For ex- 
ample, the limit on the total right-handed muon coupling 



J|.9^J 2 + |3LJ 2 + |3^| 2 + 3| 5 LJ : 



is reduced by a factor of six from the pre-TWIST value 
to <8.2xl0~ 4 (90% C.L.). 

Left-right symmetric models extend the SM with a 
right-handed W [1]. In the generalized (or non-manifest) 
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